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We consider a strongly interacting quasi-one dimensional Bose gas on a tight ring trap subjected 
to a localized barrier potential. We explore the possibility to form a macroscopic superposition of a 
rotating and a nonrotating state under nonequilibrium conditions, achieved by a sudden quench of 
the barrier velocity. Using an exact solution for the dynamical evolution in the impenetrable-boson 
(Tonks-Girardeau) limit, we find an expression for the many-body wavefunction corresponding to a 
superposition state. The superposition is formed when the barrier velocity is tuned close to multiples 
of integer or half-integer number of Coriolis flux quanta. As a consequence of the strong interactions, 
we find that (i) the state of the system can be mapped onto a macroscopic superposition of two Fermi 
spheres, rather than two macroscopically occupied single-particle states as in a weakly interacting 
gas, and (ii) the barrier velocity should be larger than the sound velocity to better discriminate the 
two components of the superposition. 
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I. INTRODUCTION 

Macroscopic superpositions are at the heart of quan- 
tum information devices as they realize quantum two- 
level systems (Qubits). Qubits can either be of single par- 
ticle nature (single atom, single spin), or realized using 
collective degrees of freedom whose low-energy spectrum 
reduces to two discrete states (eg collective internal state 
transitions for Rydberg atoms and current states in 
superconducting SQUIDs Q). Collective-mode superpo- 
sitions are expected to be better protected against some 
forms of decoherence such as particle losses, since quan- 
tum correlations are spread over several single particle 
modes and the loss of one particle does not imply the 
destruction of the collective mode, hence allowing the 
superposition to survive. On the other hand, multimode 
superpositions imply the use of several different single 
particle states, and hence typically have a limited degree 
of entanglement. 

Ultracold atomic gases are interesting candidates for 
the realization of macroscopic superpositions due to their 
high purity and tunability. Most of the current proposals 
are based on two-mode Bose-Josephson junctions [^-Q. 
Experimental advances in the realization of ring traps 
[6l-[l2| make it realistic to consider other macroscopic 
superpositions, eg the (collective-mode) superposition of 
superflow states carrying different values of angular mo- 
mentum [T3l-[l6|. where the coupling between angular- 
momentum states is provided by a localized barrier which 
breaks translational invariance; an artificial gauge field 
(or rotation) 0] gives rise to tunability equivalent to 
magnetic flux in a SQUID. As a consequence of the ring 
periodicity, the energy levels of the many-particle system 
as a function of the flux $ associated with the artifi- 
cial gauge field are periodic with period $o = 2TTh/m, m 
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being the atomic mass. In an adiabatic protocol, equally- 
weighted superposition states can be realized by tuning 
the fiux near a half-integer value of the ratio ^/^o, in 
correspondence to an avoided level-crossing [l6| . 

In the case of a quasi-lD tightly-confining ring trap 
the possibility of creating macroscopic superpositions has 
been considered in detail, both for lattice and continuum 
models [l^. With respect to an ideal Bose gas, weak 
repulsive interactions induce small energy-level splittings 
which are harmful to the superpositions. Quite inter- 
estingly, in the strongly-interacting limit of impenetrable 
bosons (or Tonks-Girardeau limit) this drawback is over- 
come. Due to its fermionized character [l^, the Tonks- 
Girardeau (TG) gas displays the same energy splittings 
as for noninteracting bosons [l^. Moreover, due to 
its impenetrability, two- and three-body losses are sup- 
pressed in a TG gas [l^ , thus eliminating one of the main 
sources of decoherence in ultracold gases. The TG gas 
is therefore a very promising candidate for the realiza- 
tion of macroscopic superpositions of current states, as 
can be obtained, eg, by setting into motion a localized 
barrier potential. In the TG limit, creation of station- 
ary superposition states with velocity and 2T:h/mL, L 
being the ring circumference can be obtained by an adi- 
abatic switching on of the barrier to a velocity -Kh/mL 
corresponding to half a Coriolis fiux quantum [16| . Such 
superpositions have maximal useful correlations for in- 
terferometry [2^ [23|, with applications to ultra-precise 
atomic gyroscopes. 

We focus in this work on a sudden switch on of the 
barrier motion to velocity v multiple of irh/mL. We 
find that it gives rise to Rabi-like oscillations between 
states with velocity components and 2v, similarly to 
what has been observed for a superlattice ring 2l|. At 
specific times, we find an equally- weighted macroscopic 
superposition of two multiparticle states with velocity 
components and 2v. In view of the multimode as- 
pects of the strongly correlated TG state, several ques- 
tions are open regarding such a novel superposition, in 
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particular on its nature and on its degree of entangle- 
ment. For example, similarly to a Fermi gas, impenetra- 
ble bosons cannot occupy the same single particle state, 
hence, different from Ref.[2l!|, the macroscopic superpo- 
sition is not expected to be close to the usual "NOON" 
state, \NOON) cx [(foj)^ -I- {b\^^)^]\vac) where all the N 
atoms occupy the q — oi the qq ~ 2mv/h single parti- 
cle level, fej being the creation operator of a boson with 
momentum q. 

These questions can be addressed by our fully micro- 
scopic, analytical solution for the dynamical evolution by 
mapping onto a Fermi gas [13] • We obtain the exact ex- 
pression for the wavefunction of the superposition state, 
which schematically reads l^*) oc An^kF<k<kF[c\ + 
•^I+go ^"''-fe+go]l^°^'')' '^here is the creation operator of 
a fermion with momentum k, A is the mapping function 
from fermions to bosons and kp — nN/L is the Fermi 
wavevector of the mapped Fermi gas. Such an entan- 
gled, correlated many-body superposition is only acces- 
sible through an out-of-equilibrium drive. Furthermore, 
our microscopic approach allows to obtain the time of 
formation of the superposition, and to set constraints on 
the excitation process. Finally, it allows to simulate the 
time-of-flight signal which is the standard probe used in 
experiments. 



II. EXACT DYNAMICAL SOLUTION FOR 
NONADIABATIC STIRRING 

We consider N impenetrable bosons on a ring of cir- 
cumference L at zero temperature, subjected to the stir- 
ring delta-barrier potential U{x,t) — Uo5{x — vt). The 
Hamiltonian is 



N 

E 



2m dx^ 



U{x,,t) 



^gSixj -Xi) (1) 



where in particular the Tonks-Girardeau regime corre- 
sponds to the limit g ^ oo and the interaction term 
can be replaced by the condition that the many-body 
wavefunction vanishes at contact between each pair of 
bosons, "^^Bi-'-Xj = xg...) = 0. As we want to describe 
the ring geometry, we impose periodic boundary condi- 
tions ie \E'b(...Xj...) = ^B{---Xj + L...) for any j = l...iV. 
The exact solution for the many-body wavefunction is ob- 
tained by the time-dependent Bose- Fermi mapping [23 |. 



^Bixi,...XN,t) ^ Ail/VNl)det[Mxm,t)] (2) 

where A = Ilj^iSgii{xj — xg) is the mapping function, 
with sgn(a;) = 1 for a; > and sgn(a;) = — 1 for a; < 0. 
The orbitals ij^i (cc, t) are obtained from the solution of the 
time-dependent Schroedinger equation 





FIG. 1. (Color online) Single-particle excitation spectrum 
kjL/ii versus stirring momentum qL/n. The color code de- 
notes the branches of given angular momentum in the zero 
barrier limit. The inset shows the short-time evolution (time 
in units of to ~ mL"^ /-kK) of the integrated particle current 
(in units of -kHN /mL^). 



We choose as initial condition for the TG gas its ground 
state in the presence of a nonmoving barrier, ie a Fermi 
sphere for the mapped Fermi gas built with the orbitals 
ij^i (a;, 0) = (j^i^^ {x) , eigenvectors of the Schroedinger equa- 
tion for a nonmoving barrier with eigenvalues Ef'^ = 

h^kf'^ /2to, for / = l...A^. The barrier is then suddenly 
set into motion at time t = 0^. The use of two uni- 
tary transformations Ui — ^ = e™""^/^ maps 
the problem onto a stationary one with twisted boundary 
conditions, enabling to express the solution of Eq.([3]) as 

^i{x,t) - e^«-e-^«'*/2™^c,7e-^^^V,(^-fi), (4) 
j 

where we have defined the quasimomentum q = mv/h. 
The orbitals 4>j{x) are the solutions of Ej(f>j(x) — 
{—{h?/2m)d^ + Uo6{x))(j)j{x) with twisted boundary con- 



ditions <f>j{x + L) = e ^'^^(l>j{x), and read 



1 j e*«'5'(e*('°^(^+'5'^ + AjC^'^'^^^^^^^) 



in [-f,0) 
in [0,|], 
(5) 



with normalizations Mj 
amplitudes Aj — sin[(fcj 



L(l 



^1 



sin{kjL) - 



q)L/2]/ sin[{kj - q)L/2], and 
wavevectors kj given by the solution of the transcen- 
dental equation kj — X sm{kj L) / {cos{qL) — cos(fcjL)), 
which determines the energy eigenvalues Ej = h^k'j /2m, 
with A = mUo/h'^. Information about the initial condi- 
tion enters Eq. (|4]) through the time independent overlaps 



Cjl 



III. 



STATE OF THE SYSTEM UNDER THE 
STIRRING DRIVE 



ihdtipiix,t) = 



-^dl + Uo6{x-vt)]Mx,t)- (3) 
2m ' 



We first consider a velocity close to a special value 
V = hmr/m.L with integer n, which corresponds to 
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half integer values of the Coriolis flux ratio $/$0i with 
^ = vL. We focus on the small barrier limit XL < 1. 
As is seen in FiglU this choice of v corresponds to sev- 
eral avoided level crossings of single particle states. As 
a consequence of the sudden quench of the barrier ve- 
locity, we find that the TG gas, which initially occupies 
the zero-momentum Fermi sphere of the mapped Fermi 
system, oscillates between two A^-particle Fermi spheres, 
one centered at fc = and the other at fc = 2q, realiz- 
ing at half oscillation an equally weighted macroscopic 
superposition of the two Fermi spheres. 

The derivation reads as follows. According to the ex- 
pression for the overlaps cji , the states excited under the 
effect of the stirring drive are fixed by quasimomentum 
conservation kj — k^^^ ± q for v < vp, or kj = ifej'"'' + q 
for V > vp = Ntm/mL. In detail, taking for simplicity 
V > vp which will turn out to be the most favourable 
situation, we find that to leading order in \L only four 
states j are coupled to each single-particle level I of the 
initial-state Fermi sphere, with coefficients \cji\ = 1/2 for 
j = n±2Int[//2] and j = n + l±2Int[l/2], Int[..] denoting 
the integer part; with the exception of the lowest state 
1 = 1, where \cji\ = 1/V2 for j = n and j = n + 1. 
For each level kj we know its momentum (hence angular 
momentum) components from the analysis of the zero- 
barrier limit, where a true level crossing occurs of two 
states of well defined angular momentum. For example, 
for V = ATrH/mL and N = 3, the level Z = 1 is coupled 
by the stirring barrier to the states with j = 4 and 5 
which are both an equal-weight superposition of states 
with momentum k = and k = 8n/L. Similarly, the 
levels I = 2 and 3 yield an equal-weight superposition 
of states with k = 2tt/L and Qii/L for the lowest-energy 
doublet and of k = —2tt/L and IOtt/L for the highest- 
energy doublet. Summing up all the contributions, we 
find that each momentum state has the same occupa- 
tion, and the momentum occupation distribution is a su- 
perposition of the two Fermi spheres {— 27r/L, 0, 27r/L} 
and {67r/L, Stt/L, IOtt/L}. A similar reasoning holds for 
arbitrary barrier velocities and particle numbers (cho- 
sen odd to ensure proper boundary conditions on the 
mapped Bose gas), leading to the occupation of two 
Fermi spheres centered at and 2q. From Eqs.([2]) and 
, the many-body wavefunction of the superposition is 
finally obtained by mapping onto a Fermi gas where each 
atom occupies a different superposition of a few (typically 
four, in the weak barrier limit) single particle orbitals. 

The detailed dynamics of the system is also sim- 
ply described, according to the values of the overlaps 
Cji, in terms of the occupation of a few momentum 
states for each single-particle state I. Using Eq. 
the time dependence eg of the particle current density 
j{x,t) = {h/m)lin^i'ipi{x,t)dxi'i{x,t) is fixed by the 
time-evolution factors e~^^^^~^i'^*/^ . The typical time 
scale is fixed by the energy level splitting associated 
to the highest two levels occupied through the stirring 
drive. The short-time behaviour of the integrated cur- 
rent I{t) = (1/i) / dxj{x,t) is illustrated in the inset of 



FiglU the time evolution is not purely sinusoidal due to 
the multimode nature of the superposition. 

Consider now the off- resonant case h{n — Vji: /mL < 
V < hmr/mL. In this case a weak barrier does not trans- 
fer angular momentum to the gas and the momentum 
occupation distribution is a single Fermi sphere centered 
at fc = 0. This is readily derived by inspecting the over- 
laps Cji, yielding \cj^i\ = 1 for j = n, and \cji\ = l/-\/2 
for / > 1 with j = n ± 2Int[//2]. 

We note that the state of the system found under a 
sudden switching on of the barrier velocity is very differ- 
ent from the one obtained by an adiabatic turning on, 
where the lowest N single particle energy levels are pop- 
ulated and the momentum occupation distribution is ei- 
ther a single Fermi sphere centered at k = q for even 
n, or a superposition of two Fermi spheres centered at 
k = q — hn / L and k = q + hn / L for odd n. 

As a partial summary, we find that the creation of 
macroscopic superpositions is efficient only in the vicinity 
of u ~ hmr/mL. Fluctuations on the barrier velocity 
might degrade the quality of the superposition, leading 
in particular to a different weight for the two components 
of the superposition. The velocity window useful for the 
excitation depends on the details of the barrier, eg on 
the barrier height. We find that an increase of the barrier 
height improves the robustness of the superposition state 
with respect to barrier velocity fluctuations. 

IV. ONE-DIMENSIONAL MOMENTUM 
DISTRIBUTION AND WIGNER FUNCTION 

We illustrate the formation of the macro- 
scopic superposition by following the dynami- 
cal evolution of the ID momentum distribution, 
n{k,t) = J dx / dye**''(^^^'pi(a;, y, t). It is defined in 
terms of the (time-dependent) one-body density matrix 
Pi ( x, t), which can be efficiently calculated according 

to [ig 

N 

pi{x,y,t)= ^ i}*{x,t)Ai^i,{x,y,t)i}v{y,t), (6) 

with Ai,r(x,y,t) = [Aet P]{P~^)Ji, and Pi,i,{x,y,t) = 
Si^i' ~ 2 dx''ipl{x',t)'ipi'{x',t). As shown in FigH dur- 
ing the time evolution induced by a stirring velocity close 
to an integer multiple of hn /mL the momentum distribu- 
tion evolves from a single peak at fc = to a single peak 
at fc = 2g, displaying at intermediate times a double peak 
structure, reflecting the superposition of the two Fermi 
spheres of the mapped Fermi gas. The peaks in the mo- 
mentum distribution, associated with the bosonic nature 
of the gas, allow to well identify the two components. 
However, since the width of the TG momentum distri- 
bution is the same as the fermionic one, in order to bet- 
ter resolve the superposition one needs stirring velocities 
larger than twice the Fermi velocity ie the sound velocity 
of the TG gas (see the inset in Figl2] for an exemple of 
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FIG. 2. (Color online) Time-dependent momentum distribu- 
tion (left panels, in units of L), as a function of the wavevec- 
tor kL/n, and corresponding Wigner function (right pan- 
els, dimensionless) as a function of X/L and kL/n at times 
t/to = 0, 2.46 and 4.92 for N=9 particles and stirring veloc- 
ity V = 14ft7r/mL. The inset of the third panel shows the 
momentum distribution at the time of an equally-weighted 
superposition for stirring velocity v = hTr/mL 



momentum distribution for a superposition state in the 
case of a small velocity v = nh/mL). 

The observation of a double peak structure in the mo- 
mentum distribution does not necessarily imply the ex- 
istence of quantum correlations between the two Fermi 
spheres [1^ . In order to evidence the nonclassical nature 
of the macroscopic superposition we compute its Wigner 
function, 

fw{X,k,t)^ [ dre'''-pi{X + r/2,X-r/2,t). (7) 



At the time of the equally-weighted superposition, the 
Wigner function displays some negative regions, see again 
FigUl This illustrates the quantum correlations between 
the two Fermi spheres, which could be quantified eg fol- 
lowing [131. 



V. OBSERVABILITY IN ULTRACOLD ATOMIC 
GASES 

Current experimental detection techniques are based 
on time-of-flight (TOF) images, obtained by releasing the 
confining potential. Neglecting the effects of interactions 
after the release from the trap, we describe the spatial 
distribution of the atomic cloud after expansion as the 
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FIG. 3. (Color online) TOF images in the k^ ky plane (in units 
of tt/L) for A'^ = 9 bosons stirred at velocity v — lAirh/mL, 
at times t/to = 0, 2.46, 4.92, from top to bottom. 



momentum distribution of the gas before expansion. In 
the case of the ring we have 



where the expression in cylindrical coordinates of the 
3D one-body density matrix on a tight ring trap of 
radius R is p™»(r, 6*, z; r', 6*', z'; t) = S{r - R)b{r' - 
R)6{z)d{z')pi{Re,R0',t). Figure [3] shows the TOF im- 
ages corresponding to various stirring times, illustrating 
the transition between a zero-current state at initial time 
to a state of angular momentum C/N — 2mvR. The ini- 
tial peak at k = deforms spirally and finally tends 
to a ring, the latter in agreement with the predictions 
of [13] for a state with well-defined current. Note that 
the TOF image of the equal-weight macroscopic super- 
position, represented in the second panel in FigjSl is not 
simply obtained as a combination of the TOF images of 
well-defined current states (first and last panel in Figl3]) , 
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due to interference between the zero-current state and 
the state at velocity 2v. 

VI. SUMMARY AND PERSPECTIVES 

We have studied the nonadiabatic excitation of 
strongly interacting bosons on a ring by a sudden set 
into motion of a localized barrier potential at velocity v. 
If the velocity is suitably chosen, the state of the system 
oscillates between a zero-current state and a state of ve- 
locity 2v, displaying a superposition of the two states at 
intermediate times. Due to the strong interactions, the 
nature of the superposition state is very different from 
the one of a weakly interacting gas; we find a superpo- 
sition of two Fermi spheres rather the superposition of 
two macroscopically occupied single-particle levels. The 
superposition of current states is evidenced by a double- 
peak structure in the momentum distribution. Due to 
the underlying multimode nature of the state, the two 
peaks are resolved at best for stirring velocities larger 



than twice the sound velocity. We have also verified the 
nonclassical nature of the superposition by the study of 
the Wigner function. Our results confirm the TG gas 
as a promising candidate for applications to quantum- 
limited metrology and seem accessible to state of the art 
experiments on ultracold gases. In perspective, it will be 
important to develop detecting techniques capable to es- 
timate the relative weight of the two components of the 
superposition. A measure of their coherence could be in- 
ferred by extending the full-counting method proposed 
in 0. 
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